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ABSTRACT 

The temperature structure of Be star circumstellar disks at the sub-solar 
metallicity appropriate to the Small Magellanic Cloud (SMC) is investigated. It 
is found that for central stars of the same spectral type, Be star disks in the 
SMC are systematically hotter by several thousand degrees compared to Milky 
Way (MW) disks with the same density structure. For early spectral types (BOe 
- B3e), this results in systematically smaller Ha equivalent widths for Be stars 
in the SMC. The implication of this result on Be star frequency comparisons 
between MW and SMC clusters is shown to be a 5 - 10% lowering of the detection 
efficiency of Be stars in SMC clusters. These calculations are also compared to 
the known Ha equivalent width distributions in the MW and SMC. For the MW, 
reasonable agreement is found; however, for the SMC, the match is not as good 
and systematically larger Be star disks may be required. 

Subject headings: stars: circumstellar matter - stars: emission line, Be - galaxies: 
Magellanic Clouds 



Introduction 



Be stars are defined observationally as non-supergiant B stars that have, or h ave had in 
the p ast, emission detected in the Balmer series of hydrogen, most notably in Ha ( jSlettebak 
19881 ). The emission is thought to ar ise from a thin, equatorial disk of gas surrounding the 
central star (IPorter fc Riviniusl 120031 ). This picture is consistent with many other properties 
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of Be stars, such as an excess o 



polarization (IMcLean fc Brown 



that directly resolve the disks ( 



1978 



infrared radiation flCote fc Waters! 119871 ). linear continuum 



Poeckert et al 



Quirrenback et al 



1997 



1979). and interferometric observations 



Tvcner et all 120051 ) 



In the Milky Way about 17% of all non-supergiant B stars are Be stars , although the 
fract ion varies widely with spectral type, reaching a maximum of ~34% at Bl (IZorec fc Briot 



19971 ). Why some B stars, and not others, become Be stars is currently unc lear. Rapid 
rotation of the central B star seems to play a key role (IPorter fc Riviniusl 120031 ). but details 
are still lacking as a definitive determination of the actual rotation r ates of the centra l stars 



Cranmer 


2005; 


Fremat et al. 


2005) 



Also unclear is the exact mec hanism(s) that creates Be s t ar disks. Keplerian rot ation 



now established for Be star disks (jHummel &: Vrancken 



gests viscous disks as the correct physical model (ILee et al 



2000 



mate rial is feed into the inner edge of the disk is unknown ( 



201ll ). It is possible pulsation plays a criti cal role (ICranmerl 120091). and there is tantalizing 



observation evidence supporting this view (IRivinius et al. 



Oudmaiier et al\ 120081) . sug- 



1991 



Cranmer 



1998 



of binarity in the Be phenomena is unclear ( Porter fc Riviniusl 2003h. although the majorit y 



Porterill999h. but how 



200 



9: 



Krticka et al. 



Huat et all 120091 ) . The role 



of Be stars do not seem to be the result of binary evolution (IVan Bever fc Vanbeverenlll997l ). 
In this current work, we consider only single, isolated Be stars. 

Key to understanding the factors that produce the Be phenomenon are cluster stud- 
ies which allow ages to be assigned to the individual Be stars. Then any trends with 



nomenon occurs primarily in the 


ater half of the main sequence ( 


Fabreeat & Torre ion 


2000; 


Fabreeat & Gutierrez-Soto 


2005: 


McSwain & Gies 


2005: 


Martavan et al. 


2006) or thr 


ough- 
s find 


out the entire main sequence ( 


Keller et al. 


1999; 


Martavan et al. 


2007a 


). Other studie 



Be stars both in very young clusters, < 10 Myr, and an enhanc ement in older clusters, 
20 — 30 Myr (IWisniewski fc Bjorkmanl 120061 ; iMathew et al.l 120081 ). suggesting two effects 
at work: some Be stars are born as rapid rotators while others become rapid rotators 
during main sequenc e evolution due to the internal redistribution of angular momentum 
jEkstrom et al.ll2008h . 



Be stars are also detectable in open clusters in the Large and Small Magellanic Clouds 
(LMC and SMC), and this allows investigation of the role of metallicity in the Be phenomena. 
The solar metalli city is z = 0.02 w hile the metallicity of the LMC is about 0.007 and 
the SMC, 0.002 (IMaeder et al.lll999l ). The trend of the fraction of Be stars in a cluster 
with metallicity seems clear from many studies: the fraction of Be stars t e nds to rise with 
decreasing metallicity ( IMaeder et al.l Il999l ; IWisniewski fc Bjorkmanl 120061 ; iMartayan et al. 
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20061 . l2007ai |2010| ). Such a trend may have a natural explanation in the depen dence of 
mass-loss due to stellar winds on the metallicity (IVink et al.ll200ll ; iPenny et al.ll2004l ); a lower 
metallicity results in weaker stellar winds and hence st ars in low-metallici ty environments 
rotate more quickly, a fact well established for the SMC ( IHunter et al.ll2008l ). Thus the more 
quickly rotating B star population leads to a higher fraction of Be stars. 

Nevertheless, all of these comparative studies between the MW, LMC, and SMC proceed 
by counting candi date Be stars in the total cluster population, selected either photometrically 
(see, for example, iMcSwain fe Gies" 2005 ) or spectroscopically with very low-resolution (see 
Martayan et al. 201oh . Of course not all Be stars in a population can be found this way and 
detection efficiencies generally point to selection criteria of a Ha emi ssion equivalent widtlj j] 
in excess of 5 — 10 A as required to identify a candidate Be star ( Martayan et al. 201ol ). 
Assumed in all of these studies is that individual Be stars are more or less the same in the 
MW, LMC, and SMC so that the comparative statistics are not affected by the detection 
efficiency However, there are indications that, at some level, this may not be true. For 
example, the equivalent width distribution of Ha differs b etween the MW and S MC, with 
the SMC distribution peaking at higher equivalent widths (IMartayan et al.ll2007bl ). 



In this work, we will illustrate a fundamental difference between Be star disks in the 
MW and the SMC: for the same central star spectral type and disk density distribution, disks 
in the SMC are much hotter than in the MW. We will quantify this effect over the early Be 
star spectral types (those for which LMC and SMC samples are reasonably complete), and 
we will examine how counts of the Be star populations based on the Ha equivalent width in 
different metallicity systems might be affected by this systematic temperature difference. 



2. Calculations 

The thermal str ucture of the Be star disks was computed with the bedisk code of 
Sigut fe Jonesl (120071 ). This code enforces radiative equilibrium in a photoionized circum- 
stellar disk including the heating and cooling processes for the 9 most abundant elements (H, 
He, CNO, Mg, Si, Ca, & Fe) over several ionization stages. Details of the atomic mo d els an d 
atomic data, as well as an overview of the bedisk code, are given in lSigut fe Jonesl (120071 ). 



The density structure of the disk was assumed to be of the fixed form 



p(R,Z) 



(1) 



We adopt the convention that a positive equivalent width denotes emission. 
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where R and Z are the cylindrical co-ordinates for the assumed axisymmetric disk, and R* 
is the radius of the central B star. The quantities n and p Q are adjustable parameters that 
fix the density structure of the disk. The function H, defined below by Eq. [21 sets the 
vertical (or Z) scale height of the disk. This simple de nsity model has been very successful 



in interpreting a wide range of Be star observations (IGies et al.l 120071 ; iTycner et al.l 12008 



Jones et al.l 120081 ). Thus the thermal models of this work are identified first by the spectral 



type of the central B star (which sets the photoionizing radiation field) and then by the disk 
density parameters p Q and n. 

The fundamental parameters assumed for the central B stars are given in Table [TJ 
The mass and ra dius as a f unction of spectral type, as well as the MW T e g calibration, 
are adopted from ICoxi (120001 ) . Note that we have assumed the same mass and radius as a 
function of spectral type for the MW and SMC. Stellar evolution calculations have shown 
that B stars are likely sm aller in low metallicity systems such as the SMC. For example, 



Maeder fc Meynetl (120011 ) show that a 20 M star (roughly a B0 star) decreases in radius 
by about 20% as the metallicity is lowered from Z = 0.02 to Z = 0.004. However, we 
have ignored these differences to ensure both the SMC and MW disk models have the same 
density structure. The spectral-type - T eff calibration for the SMC is discussed below. 

The vertical (or Z) dependence of Eq. [1] contains the function H which is defined as 



2R 3 kT d 



1/2 



H = (2) 
[GM^m n \ ' 1 ; 

Here M* is the mass of the central B star, /i is the mean-molecular weight of the disk gas, 
and Td is an assumed isothermal disk temperature. Eqs. [JJ and [2] follow from the assumption 
that the disk is in vertical hydrostatic equilibrium set by the Z component of the stellar 
gravitational acceleration and an assumed isothermal disk temperature, Td. Note that T d is 
only used in Eq. [H to fix the vertical density structure of the disk; the actual temperatures in 
the disk, T(R, Z), are found by enforcing radiative equilibrium. 



Si gut et al.l (120091 ) consider consistent models in which the vertical structure of the 



disk is found by integrating the equation of hydrostatic equilibrium in a manner consistent 
with the radiative equilibrium disk temperatures; this treatment eliminates the need for the 
parameter T^. However, the disk density distribution is now dependent on the temperatures 
in the disk. As the desire of the current work is to compare the temperature of MW and 
SMC disks of identical density structure, we have used Eqs. [1] and [2] with the same Td for 
each spectral type (namely 0.6T e g W where T e ^ w is the Milky Way effective temperature at a 
given spectral type)E| We do not expect this assumption to significantly affect our differential 



2 Their remains a difference between the MW and SMC density structures due to the different mean- 
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comparison. 



The disk is assumed to be in Keplerian rotation and hence rotationally supported in 
the R direction. Gravitational darkening due to the potential rapid rotation of the central 
B star is not inclu ded, although this process can affect the thermal structure of the disk 
jMcGill et alJboilh . 



The main energy input into the disk is assumed t o be photoionizing radiation from 
the central sta r . For this paper, the older iKuruczl ( 119931 ) LTE photoionizi ng fluxes used by 
Si gut fc Jonesl ( 120071 ) were replaced with the newer non-LTE calculations of lLanz &: Hubeny 
( 120071 ) . Two grids from these calculations were used, the solar metallicity grid (pOO) to 
compute the MW models, and the 1/10 solar abundance grid (mlO) to compute the SMC 
models. 

The thermal disk models are thus described by three parameters: the spectral type of 
the central B star (which is assumed to fix the stellar mass, radius and effective temperature), 
and the two parameters in Eq. [T] that fix the density of the disk, p a , the base density (in 
g cm -3 ) and the power-law index n. For this set of calculations, spectral typ es B0, B0. 5, Bl, 



B1.5, B 2, and B3 were considered as these are the spectral types for which iMartayan et al. 



( j2007af ) consider their SMC cluster samples complete. Power-law indexes n 
and 11 densities ranging from 10 -13 to 10~ 10 gem -3 were considered. 



2, 3, and 4 



Two sets of Be star disk models were computed, one appropriate to the MW and one 
appropriate to the SMC. Three basic changes distinguish between the MW and SMC models: 

[1] The abundance table used as input to bedisk was changed to reflect the composition 
of the gas. Th e adopted abundances for the MW and SMC are given in Table The solar 



abundances of lAsplund et al.l ( 120051 ) were assumed for the Milky Way, and the abundance 



2007 



2009). The adopted SMC abundances agree very well with the recent compilation of 



Rolleston et al. 


2003; 


Hunter et al. 


2005; 


Evans et al. 


2005: 



Evansl (120091 ). The general abundance pattern of the SMC is A[Z/H] ~ —0.7 ( Mokiem et al. 



20071 ) but we keep the individual abundances for each element. 



[2] The input photoionizing radiation field was taken to be either the lLanz fc Hubeny 



( 120071 ) pOO for MW runs or mlO for the SMC runs. Test calculations at the solar abundances 
of Table H] show that very similar temperatures result using either the pOO (solar) or mlO 
(1/10 solar) photoionizing radiation fields. Hence, the small mismatch bet ween the overall 
A[Z/H] ~ —1.0 abundances used to compute the photoionizing radiation by lLanz &; Hubeny 



molecular weights of the gas, but this difference is extremely small. 
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( 120071 ) and actual SMC abundances used here is likely of little consequence to the results. 



[3] It is well established that the spectral-type - T eS ca libration differs bet ween the 



MW and the SMC, with the SMC early-B stars being hotter (bundle et al]l2007h. Thisis 



shown in Figure [TJ Here the M ilky Way ca libration of B dwarfs by iTrundle et al.l (120071 ) 
agrees well with the T e g scale of ICoxl (120001 ) (adopted in this work) whereas the LMC and 
SMC calibrations are systematically hotter by a few tho usand degrees. Our a dopted T eS 
scales are given in Tabled) Note that we do not adopt th e ITrundle et al.l (120071 ) SMC T ef j's 
dir ectly but increase ou r adopted MW T eS calibration of ICoxl J2OO0F by the amount found 
by ITrundle et al.l (120071). Howe ver, the differences between our SMC T e ff's and that found 
directly by ITrundle et al.l ( 120071 ) are small. 



Of these three differences between MW and SMC models, [1] and [3] have the largest 
affect on the disk temperatures. 

One caveat to our analysis is that we have ignored gravitational darkening of the central 
B stars due to rotation. As a class, Be stars a re well established to be rapid rotators 
(jPorterl Il996l : lYudinll200ll ; iPorter fc Riviniusl 120031 ). although how close to critical rotation 
(in which the equator ial velocity of the star equals the orbital speed) is still contentious. 
Townsend et al.l (120041 ) demonstrated that rotation rates above 80% of critical rotation (or 
Wfrac > 0.8) may be difficult to detect because gravitational darkening lead s to a significant 
cooling of the star's most rapidly rotating equatorial regions. Nevertheless, ICranmerl (120051 ) 
used a statistical analysis including the effects of gravitational darkening to demonstrate that 
the distribution of rotation speeds for early Be stars (in the MW) is fairly uniform between 
60 and 100 % of critical. Hence for the MW, a representative value of 80%, or t>f rac ps 0.8, is 
reasonable. iMartayan et al.l ( j2007al ) demonstr ate that for t h e SM C field NGC 330, Be stars 
rotate on average at t>f rac = 0.87. Recently, iMcGill et al.l (120111 ) have demonstrated that 
Be star disks become systematically cooler as a result of gravitational darkening. Thus any 
potential difference in mean rotation rates, t>f ra c, between the MW and SMC could result in 
a system atic decrease in the temperature of SMC Be star disks. Fortunately, we can use the 
results of lMcGill et al.l ( 1201 ll ) to estimate the expected effect. Taking t>f rac = 0.8 for the Milky 
Way and t>f rac = 0.9 f or the SMC (roundi ng the observed value up for the maximal effect), 
Figures 4, 5 and 6 in IMcGill et al.l ( 1201 ll ) give a reduction in the density-weighted, global 
disk temperature of ~300 K for an increase in Vf rac from 0.8 to 0.9. As will be demonstrated 
in the subsequent sections, this temperature decrease is much smaller (by at least a factor of 
five) than the temperature increase predicted for SMC disks due to effects [1] and [3] above. 
As a result, the neglect of gravitational darkening is not a significant source of uncertainty. 
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Fig. 1. — Comparison of the spectral type -T P ff calibr ation between the MW and the 
SMC/LMC. Shown a re the MW calibration of Icoxl (hoOflh and the LMC/SMC calibra tion of 



Trundle et al.l (120071 ) . The x-axis is the s- index of Ide Jager fc Nieuwenhuijzenl ( 119871 ) . The 



corresponding spectral types are indicated at the top of the figure. 
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3. Results: Temperature Structure 

In this section, we consider how disks of the same, fixed density differ in their thermal 
structure due to the metallicity of the gas. In order to compare a large number of models of 
differing spectral types and density parameters, we will construct a single, global measure 
of the disk temperature, namely the density- weighted, average disk temperature, defined as 

<T> disk =-^- / T(R,Z)p(R,Z)dV . (3) 
^ disk J disk 

Figure [2] shows this density- weighted disk temperature as a function of the central B star's 
T e ff for both the MW and SMC for all density models with a radial power-law index of 
n = 3 in Eq. [TJ Figure [3] shows a similar plot but for the slower radial drop-off of n = 2. 
In both figures, the global disk temperature is shown both as an absolute temperature and 
as a fraction of the stellar T e g (which differs between the MW and SMC). For each T e g, 
temperatures and temperature ratios are plotted for each of the 11 disk base densities, 
p , considered. Immediately evident from the figures is the significantly lower global disk 
temperatures for high values of p„. Dense disks develop a cool equatorial zone near the star 



(IMillar fc Marlboroughlll998l : ISigut et al.ll2009l ) and this can significantly lower the density- 



weighted average disk temperature. This effect is particularly noticeable in Figure [3] which 
includes the densest disks considered in this work (n = 2 and p Q = 10~ 10 gem -3 ). 



For n = 3, the MW temperature ratios are reasonable well fit by <T> disk /T cS « 0.5qj 
wheres the SMC disks are better fit by < T >di S k /T e g ~ 0.61. Thus the SMC disks are 
even hotter than expected by the increased effective temperatures of the SMC B stars. This 
reflects the lower cooling rates in the SMC disks due to the lower metallicity of the gas. Note 
that the densest disks (the largest symbols) are not well-fit by the general median trend and 
the fit is particularly poor for the densest disks with n = 2 (Figure [3]); the cool equatorial 
zone significantly lowers the density-weighted average temperatures. 



This entire situation is summarized in Figure H] which shows the disk temperature ratios, 
averaged over p Q , as a function of T e g for the three power-law indexes considered, n = 2, 3, 
and 4. There is little difference between the n = 3 and 4 models whereas the n = 2 models 
are significantly cooler due to their larger densities. Note the median ratio over all T e g 
considered is always higher in the SMC. The numerical values for the median temperature 
ratios are given in Table [31 



Sigut et al.l (|2009l ) suggests a value of 0.6 as most appropriate for Be star disks. However, the models in 



that work do not have a fixed vertical density structure but self-consistently solve the equation of hydrostatic 
equilibrium. In addition, they are of an odd type compared to the present work: pure hydrogen/helium disks 
(i.e. zero metallicity) with a solar metallicity (pOO) photoionizing radiation field. 
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Fig. 2. — The density- weighted average disk temperature for the MW (left panels) and SMC 
(right panels). All models assumed n = 3 in Eq. [T]and the symbol sizes indicate the size of 
p Q . The median at each spectral type is indicated by cross and the straight lines are fits to 
the medians. 
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Fig. 3. — Same as Figure [2] except all models assumed n = 2 in Eq. [TJ Note the very low 
temperatures predicted by the highest density models. 
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Figure |5] shows as a contour plot of the difference between SMC and MW density- 
weighted average disk temperature as a function of the disk density parameters p a and 
n for the spectral type Bl. Differences of up to 4500 K are seen (SMC hotter) with the 
largest values occurring for intermediate disk densities p Q . The location of the maximum 
temperature difference in this plot tends to move to higher p Q for larger n. The temperature 
difference rapidly drops for the densest models considered (logp G > —10.4 and n < 2.5) as 
the density-weighted temperature in this region is dominated by the very cool equatorial 
zone close to the star; this occurs in the most optically thick portion of the disk where the 
sensitivity to the gas metallicity is least. 

Figure [6] again shows the temperature difference between the SMC and MW models but 
this time as a function of spectral type and p n for a fixed power-law inde x of n = 3. The 
spectral type is plotted as the s-factor of lde Jager fc Nieuwenhuijzenl (119871 ) (namely 1.20 for 
BO, 1.35 for BO. 5, 1.50 for Bl, 1.65 for B1.5, 1.80 for B2 and 1.95 for B3). Again the largest 
temperature differences are for the intermediate density models with the largest absolute 
differences occurring near spectral type B1.5, near the peak of the Be star distribution with 
spectral type. 



4. Results: Ha Equivalent Width 



While interesting, the previous results are of limited practical use because the thermal 
structure of a circumstellar disk is not directly observed. The main observable for Be stars 
is the equivalent width (in emission) of the Ha line (A 6562.8 A) in the spectrum. Hence 
it is of direct practical interest to compare the Ha emission line equivalent widths between 
the MW and SMC models. All current methods designed to find and measure the Be star 
fraction in a stellar cluster implement some t echnique to selec t cand idate emission-line stars 
based on a measurement of Ha. For example, iMartayan et al.l (120101 ) use very low-resolution 
slitless spectroscopy and estimate that their survey of Be stars in SMC clusters detects Be 
stars with Ha emission equivalent widths in excess of 10 A or peak intensities more than 
twice the adjacent continuum. 



To this end, we have used the beray code of lSigutl (120101 ) to compute the Ha line profile 
and equivalent width (in A) for each of the MW and SMC disk models of the previous sec- 
tion. The beray code solves the equation of radiative transfer along a series of rays through 
the star-disk system to compute both resolved monochromatic images and unresolved spec- 
tra. The Ha calculations used the n = 2 and 3 level populations of hydrogen from the 
bedisk thermal solution. Rays that terminate on the stellar surface use the Doppler shifted 
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Fig. 4. — The median (over p Q ) density-weighted disk temperatures for the MW (left panel) 
and SMC (right panel) as a function of stellar T e g . The symbols denote the power-law index 
in Eq. [T]as n = 2 (circles), n = 3 (squares), and n = 4 (diamonds). The horizontal lines are 
the best constant fits to the medians. 
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Fig. 5. — Difference (SMC-MW) in the density-weighted average disk temperature, Eq. [3J 
as a function of the disk base density, p D , and power-law index, n, for spectral type Bl. The 
contours are of equal temperature difference in degrees K. 
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Fig. 6. — Difference (SMC-MW) in the density-weighted average disk temperature as a 
function of the disk base density, p Q , and the spectral type of the star, s-factor, for power- 
law index n — 3. The contours are of equal temperature difference in degrees K. 
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photospheric LTE Ha line profiled of Barklem Sz Piskunov ( 2003 ) for the initial boundary 
condition, while rays that pass entirely through the disk assume no incident radiation. The 
formal solution to the transfer equation was performed using the Ha source function and 
the bound-free and free-free continuous opacity and emissivity sources at the wavelength of 
Ha. Electron scattering was also included using the mean intensity at each position in the 
disk as the source function for the scattering emissivity. The disk was assumed to be in pure 
Keplerian rotatio n and the freq uency- angle dependent transfer equation was solved in the 
observer's frame (IMihalaa Il978l ) . 



Calculating the Ha line profile adds two parameters to the problem, the inclination 
of the Be system to the line of sight (i = 0° indicates a pole-on star/disk, and i = 90°, an 
edge-on one) and the outer radius of the disk, Rd (quoted here in units of the stellar r adiu s )@ 
Values for these parameters were taken to be % = 15, 30, 45, 60, and 75°, and Rd = 5, 10, 
15, and 20 R*. For each spectral type, there were a total of 660 models representing all 
permutations of the parameters p a , n, Rd and i. 

The ratio of the SMC and MW Ha equivalent widths is shown in Figure [3 The trend is 
clear: Ha is generally weaker in the SMC compared to a disk of the same density structure 
in the MW. Over all models satisfying EW(Ha) > 5 A, only 366 out of 2047 models (or 18%) 
had a stronger Ha line in the SMC; the median SMC to MW equivalent width ratio is 0.89. 

For early-type Be stars, disk temperatures usually exceed 10 4 K and increased temper- 
atures tend to reduce the population of the excited states of hydrogen due to increased 
ionization. This results in a weaker Ha line. Neverthele ss, the temperatu re structure of the 
denser Be disk models is quite complex (see, for example. ISigut et al.ll2009i ) and some of these 
models can produce a stronger Ha line in the SMC. Table |4] gives the fraction of the models 
in which the SMC Ha equivalent width was larger than the MW model (for an identical 
disk) as a function of spectral type. This fraction increases strongly for later spectral types 
(more than doubling over the range considered) showing that the cooler disks around later 
spectral types have a much larger fraction of increased Ha strengths in the SMC. 



4 These profiles were computed in stellar atmospheres assuming the solar abundance. While the profile 
corresponding to the required stellar temperature (T^ w or Tj^ 10 ) was used, this still represents a small 
inconsistency. We suspect that this effect is smaller than the uncertainty introduced by the use of LTE, as 
opposed to non-LTE, profiles. Nevertheless as we are really interested only in models with strong emission 
in Ha for which the disk dominates, the adopted photospheric profile is not a major source of uncertainty. 

5 The thermal solution for the disk structure, computed by BEDISK, assumed Rd = 30i?». In beray, 
smaller choices for Rd were enforced by multiplying all level populations beyond Rd by 10~ 5 . The inner 
radius of all disks was taken to be the stellar radius. 
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Table [5] shows an additional interesting result. Over all models (spectral type and 
density structure), there is a very strong dependence of the fraction of models that predict 
a stronger Ha EW in the SMC on the viewing angle. Increased emission is essentially 
eliminated (fraction less than 10%) for models seen more "disk-on" with i > 60°. Note that 
the maximum viewing inclination in this work, 75° does not produce obvious shell-spectra 
characteristic of Be star disks seen nearly edge-onjj 

The difference in Ha equivalent widths between the MW and SMC could potentially 
affect the comparison of the fraction of Be stars in clusters. Even if both the MW and 
the SMC had the same distribution of disk density parameters (i.e. p a , n and Rd), there 
would nevertheless be a different distribution of Ha line strengths because of the different 
temperature structure of the disks. Figure M shows the SMC Ha equivalent width plotted 
against the MW value for a disk of the same density structure. If one assumes a selection 
method for candidate Be stars with the same threshold in the MW and SMC, there will 
be a population of stars with weaker emission in the SMC that fall below the detection 
limit whereas the MW counterparts would be counted. This effect is quantified in Table [6] 
which gives the missed percentage of Be stars among the SMC models as a function of 
Ha detection thresholds of between 2 and 15 A. As can be seen, this effect is small with 
the missed percentage rising to ~ 10% only for thresholds as large as 15 A. For thresholds 
between 5 and 10 A, the prediction is about 5%. Hence despite the rather large temperature 
differences between the SMC and MW disks, the systematic difference in the Ha equivalent 
widths is not a serious source of bias in Be star candidate counts. 

Finally, we compare the Ha equivalent width distributi ons predicted by our m odels 



of both the MW and SMC with the observational results of iMartayan et all (l2007bl ). In 
this comparison, we include all models producing emission (i.e. all values of n, p Q , Rd 
and i which result in EW(Ha) > 0) except those with the smallest disks, Rd = 5R*. In 
constructing the histograms, we have weighted each model based on its inclination as follows: 
for random inclinations, the probability of observing an inclination between i and i + di is 
sini di. Therefore, we have assumed that each of the five values of i considered is the centre 
of a (non-overlapping) bin and have weighted each model by the fractional area in that bin@ 



The results are shown in Figure M For the Milky Way, IMartayan et all ( l2007bl ) adopt 



6 Be shell stars are not considered in the present work. The fraction o f Be shell stars in the M W is ps 23% 
|Hanuschiklll996h . The fraction in the SMC is estimated to be w 16% (jMartavan et alj|2007bh . 

7 The bins are (in degrees) [0, 22.5), [22.5, 37.5), [37.5, 52.5), [52.5, 67.5), and [67.5, 90] which have fractional 
areas 0.0761, 0.1305, 0.1846, 0.2261, and .3827, respectively. As an example, for each 10 stars in the first 
bin, there should be 50 stars in the last bin. 
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the equivalent wi d th me asurements of lAndrillant fc Fenhrenbach I (119821 ) . lAndrillantl ( 119831 ). 
and iDachs et all (Il992l ). As noted, the only "fine-tuning" of our models is the exclusion 
of the small Rd = 5 R* disks; these disks produce a large peak for EW(Ha) < 10 A which 
is not seen in the observations. Such small disks, when they seem to occur in nature, are 
usually attributed to binary truncation (for example see the case of a Ara as discussed by 



Meilland et al.ll2009l ). Given this, the predicted MW equivalent width distribution gives a 



reasonable match to the observations. However, we do note that the sample of Milky Way 
observations finds few widths above 55 A whereas such values are represente d in our models. 



For the SMC, the models produce a poorer match to the observations of IMartayan et al. 



( 2007bl ). We suspect the large mismatch for < 10 A may be incompleteness in the observa- 
tions. However, the models under-predict the fraction of large equivalent widths (> 25 A) 
by a significant a mount and this d i fferenc e seems to be a difference between the MW and 
SMC, as noted by IMartayan et all (l2007bl ). As suggested by these authors, this and other 
observational evidence suggests that Be star disks are systematically larger in the SMC. 
Indeed, Figure [10] shows that we can produced a peak in the SMC distribution near 40 A 
by restricting the model set to only the largest disks considered, Rd = 20 R*, although the 
peak is not wide enough. We suspect that including even larger disks might improve the fit. 



5. Conclusions 

We have computed a large set of Be star disk models appropriate to the average metal- 
licities of the Milky Way and the SMC. We have shown that for a disk of identical density 
structure, SMC disks are systematically hotter than MW disks, typically by several thousand 
degrees. This difference is attributable to the higher T e g scale of massive stars in the SMC 
and the lower metallicity of the SMC disk gas, which reduces its ability to cool. For all of 
the considered models, the Ha emission equivalent width is generally smaller in the SMC 
as compared to MW disks of identical density. However, this systematic difference is not 
predicted to affect comparisons of Be star fractions between the MW and SMC. 

We also show that the Ha equivalent width distributions of our models is in reasonable 
agreement with the known MW distribution. However, a similar distributi on of underlying 



parame ters does not seem to fit the SMC. We confirm the suggestion of IMartayan et al. 



( ]2007bl ) that systematically larger Be stars disks in the SMC may explain this discrepancy. It 
would be desirable to put this interesting result on a sounder footing by directly determining 
the disk parameters for a sample of SMC Be star disks via the analysis of high-resolution 
spectra. Such a study would directly give the distribution of disk parameters (p Q , n and Rd) 
which could be compared to a similar sample in the MW and would allow a more detailed 
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Fig. 7. — Ratio of the SMC to MW Ha equivalent width as a function of the MW equivalent 
width. Only models predicting a MW equivalent width of > 5 A are plotted (see text). 
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Fig. 8. — The SMC Ha equivalent as a function of the MW Ha equivalent width. Detection 
limits of 5 A are as indicated, as are regions where disks would be detected both in the MW 
and SMC and only in the MW. The dashed line is of unit slope and the solid line is a 
quadratic, least-squares fit to the calculations. 
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Fig. 9. — Histograms of the Ha equivalent widths predicted by the MW and SMC models. 
All models predicting emission (EW > 0) are included with the exception of the models with 
very small disks (Rd = 5 R*). The fraction represented by each EW bin is rel ative to the total 



number of models satisfying EW > 0. The squares are the observations of iMartayan et al. 
(2007b|) (in the same EW bins) with y/~N error bars added. 
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Fig. 10. — A Histogram of the Ha equivalent widths predicted by the SMC models with 
20 i?*. All models predicting emission (EW > 0) are included and the frac tion represented 
by each bin is relative to the total. The squares are the observations of iMartayan et al. 
( 2007bl ) (in the same EW bins) with y/N error bars added. 
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investigation of how lower metallicity affects these circumstellar disks. 
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This work is supported by the Canadian Natural Sciences and Engineering Research Council 
(NSERC) through a Discovery Grant to TAAS. 



-24- 



Table 1. Parameters adopted the central B stars. 



Spectral Type 


Mass 


Radius 


rpMW 


j^SMC 




(M Q ) 


(Re) 


(K) 


(K) 


BO 


17.5 


7.4 


30000 


32000 


B0.5 


15.4 


6.9 


28000 


29600 


Bl 


13.2 


6.4 


25000 


28000 


B1.5 


11.0 


5.9 


23000 


26400 


B2 


9.1 


5.3 


21000 


23800 


B3 


7.6 


4.8 


19000 


22000 



Note. 



The mass and radius calibrations are 



taken from ICoxl (120001 ) and assumed to be the same 
for the M W and SMC. The MW T eS calibration is 
fro m I Cox (2000) and t he SMC calibration is based 
on 



Trundle et all (120071 ) (see text). 
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Table 2. Adopted abundances for the MW and SMC. 



Element Abundance 





MW 


SMC 


He 


10.92 


10.9 


C 


8.39 


7.34 


N 


7.78 


6.97 





8.66 


8.10 


Mg 


7.53 


6.76 


Si 


7.51 


6.77 


Fc 


7.45 


6.93 



Note. - The tabulated 
abundances are given as 
A = \og w {N x /N n ) + 12. 
The MW (solar) 
abundances are from 



Asplund et all (12005( 1 . 



Sources for the SMC are 
discussed in the text. 



Table 3. Disk temperatures ratios averaged over all spectral types, p a , and Rj, as a 
function of power-law index n. The quoted dispersion is 1 o. 



Index n 


<^>disk Meff 




MW 


SMC 


4.0 


0.57 ±0.07 


0.61 ±0.08 


3.0 


0.55 ±0.06 


0.61 ±0.08 


2.0 


0.48 ±0.10 


0.53 ±0.13 
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Table 4. Number and fraction of Be star models with a larger Ha equivalent width in the 
SMC as a function of the spectral type of the central B star. 



Spectral Type 


Ha Larger 


Ha Smaller 


Larger Fraction 


BO. 5 


38 


325 


10% 


Bl 


57 


346 


14% 


B1.5 


82 


341 


19% 


B2 


81 


346 


19% 


B3 


108 


323 


25% 



Note. — Only models predicting EW(Ha) > 5 A emission are 
included. 



Table 5. Number and fraction of Be star models with a larger Ha equivalent width in the 

SMC as a function of the viewing angle %. 



Inclination 


Ha Larger 


Ha Smaller 


Larger Fraction 


15° 


133 


332 


29% 


30° 


110 


316 


26% 


45° 


80 


347 


19% 


60° 


33 


345 


9% 


75° 


10 


341 


3% 



Note. — Only models predicting EW(Ha) > 5 A emission 
are included. 
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Table 6. Percentage of Be stars missed in the SMC as a function of Ha equivalent width 

threshold. 



Threshold EW 


Number of Detected Be Stars 


Percent Missed 


(A) 


MW 


SMC 




2.0 


2210 


2157 


2.4% 


5.0 


2047 


1939 


5.2% 


10.0 


1774 


1670 


5.9% 


15.0 


1458 


1324 


9.2% 
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